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ABSTRACT 



We present a statistical study of the environments of 63 high-redshift radio galaxies (HzRGs) between redshifts 1 < z < 5.2, using 
the 24 jim waveband of the MIPS instrument aboard the Spitzer Space Telescope. Using a counts-in-cell analysis, a statistically 
significant source overdensity is found in 1.75' radius circular cells centred on the HzRGs when compared to reference fields. We 



report an average overdensity of 5 {- Ntargets/N, 



targets I reference 



) = 2.2 + 1.2 at a flux density cut of fi 



24yum 



0.3 mJy This result implies that 

HzRGs are likely to lie in protoclusters of active and star-forming galaxies at high redshift. Over 95% of our targeted HzRGs lie in 
higher than average density fields. Further, 20 (32%) of our selected fields are found to be overdense to at least a 3cr significance, of 
which 9 are newly identified protocluster candidates. We observe a weak correlation between redshift and 24 yt/m source density, and 
discuss the populations being probed at diff'erent redshifts. In our uniformly selected sample, which was designed to cover two orders 
of magnitude in radio luminosity throughout z - 1 - 4, we find that the 24 yt/m source density does not depend on radio luminosity. 
We also compare this result with recent work describing IRAC source overdensities around the same HzRGs and find correlations 
between the results. 

Key words, clusters: observations - galaxies: evolution - galaxies: formation - infrared: galaxies - large scale structure 



1. Introduction 

Observations of high-redshift radio galaxies (HzRGs: z > 1.0 
with a rest-frame 3GHz luminosity greater than 10^^ WHz"^) 
provide evidence that such galaxies are the progenitors of the 
current day giant elliptical (gE) and cluster-dominant (cD) 
galaxies (see e.g. De Breuck et al. 2002; Seymour et al. 
2007; De Breuck et al. 2010). As such, we expect HzRGs 
to be signposts of high-density regions. Studies have shown 
overdensities of Lya emitters in HzRG environments between 
2.2 < z < 4.1 {e.g. Venemans et al. 2002, 2003; Overzier 
et al. 2008). The formation of filamentary structures out to 
high redshifts (Croft et al. 2005; Matsuda et al. 2005) and 
of Extremely Red Object (ERO) overdensities in high-redshift 
active galactic nuclei (AGN) fields have also been observed {e.g. 
Chapman et al. 2000; Hall et al. 2001). Overdensities of Ha 
emitters (Kurk et al. 2004a; Tanaka et al. 2011), Lyman Break 
Galaxies (Miley et al. 2004), selected sources (Galametz 
et al. 2009) and red sequence galaxies (Doherty et al. 2010) have 
also been reported in the vicinities of HzRGs. Large rotation 
measures within the HzRG environments (Carilli et al. 1997) are 
indicative of dense. X-ray emitting cluster atmospheres. All of 
this evidence suggests HzRGs trace protoclusters, and indicate 
that protoclusters can be in place from early epochs (when the 
Universe was only ~ 1 - 2Gyr old). These studies have been 
performed on limited samples of less than a dozen each; prior 
to this work, a systematic study of HzRG environments has not 
been undertaken. 

High-energy continuum emission from a central accretion 
disc in an AGN is believed to be re-processed by the dusty 



toroidal structure encompassing it, and re-radiated in the mid- 
IR (10 - 30yum; Elvis et al. 1994). Since all AGN are believed 
to have such a torus, they will all have a mid-IR signature, 
and though there may be orientational eflTects, all AGN are 
thought to emit strongly in the rest- frame mid-IR. However, mid- 
IR observations will not only detect AGN; starburst galaxies 
have strong polycyclic - aromatic hydrocarbon (PAH) fea- 
tures, strongest at -3.3 yum, 6.2yum, 7.7 yum, 8.6yum, 11.2//m 
and 12.7yum, attributed to vibrational modes of these complex 
molecules (Gillett et al. 1973). These attributes make mid-IR 
detection ideal for selection of star-forming galaxies at high 
redshifts as well. 

With blind targeting requiring large-area surveys, eflTorts to 
identify clusters can be both time and resource intensive. If we 
can pre- select targets that are thought to preferentially reside 
within rich environments, we can quickly identify larger num- 
bers of robust cluster candidates. Selecting HzRGs (which have 
been found to lie in overdense regions at other wavebands and at 
lower redshifts) and analysing their fields in the mid-IR allows 
us to probe the AGN and starburst populations simultaneously, 
without obvious selection eflfects. 

The companion paper to this work, Galametz et al. (2012, 
ApJ, in press, and hereafter referred to as G12), undertakes the 
analysis of 48 HzRG fields between redshift 1.2 < z < 3.0 using 
the NASA's Spitzer Space Telescope (Werner et al. 2004) Infra- 
Red Array Camera (IRAC; Fazio et al. 2004) to isolate high red- 
shift sources. Their work shows that HzRGs lie preferentially in 
medium-to-high density regimes when compared to a reference 
sample, though the colour criterion used in that work does not 
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differentiate between passive and active galaxies. 

Throughout this paper many references to the HzRG fields 
are made, they are addressed with the same IDs as the HzRGs, 
though typically we are discussing the environments and not the 
radio sources. In this paper we refer to all significantly overdense 
regions as protoclusters or protocluster candidates. Though at 
low redshifts these structures may be gravitationally bound, 
protoclusters is used throughout for the sake of consistency. A 
ACDM cosmology with Hq = 70kms"^Mpc"\ flm = 0.3 and 
Qa = 0.7 is assumed throughout. 

2. Observations 

2.1. HzRGs 

Our sample of HzRGs is taken from both flux-limited sur- 
veys such as 3C, 6CE, 7C and MRC for their unbiased radio 
properties, as well as ultrasteep radio spectra surveys (e.g. De 
Breuck et al. 2000) to cover the higher redshift (z > 2) range. 
This sample is specifically designed to mitigate the luminosity- 
redshift correlation that aff'ects flux-limited samples; the reader 
is referred to De Breuck et al. (2010, Fig. 2) for further details. 

We present observations of the first statistically significant 
sample of HzRG fields to date. The 70 target Spitzer-HzRG 
sample is described in more detail by Seymour et al. (2007) 
where targets were preferentially selected based on their amount 
of supporting data (Hubble Space Telescope and SCUBA ob- 
servations) and existing guaranteed time observations (GTOs) 
with Spitzer. With no 24 jim data (one source) and insufficient 
coverage (seven sources), our sample size is reduced to 62 
targets. The source LBDS53W069, although not meeting the 
HzRG radio luminosity selection criterion, has been included in 
this sample as it is a known radio galaxy with the relevant data 
available (Stern et al. 2006), making our HzRG sample 63 fields 
in total. 

Our HzRG sample was observed at 24yum using the 
Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 
2004). The majority of observations (59 sources) have been 
undertaken as part of Program ID's 40093 (39 sources) and 
3329 (20 sources) as part of Spitzer proposal cycles 1 and 4 
respectively. Exposure times for the former 59 sources were 
determined based on 24yum background levels as determined 
with the Leopard software. Data for three remaining sources 
were taken from much deeper GTO programs. 

2.2. Reference Fields 

The reference set against which our HzRG fields are analysed 
is that of the Spitzer Wide-area InfraRed Extragalactic Survey 
(SWIRE; Lonsdale et al. 2003). SWIRE, the widest area cryo- 
genic Spitzer legacy program, imaged 49 sq. deg. of the sky 
over six^ high-latitude extragalactic fields in all seven available 
bands^ and detected over two million objects up to z ~ 4 
(Rowan-Robinson et al. 2008). The field choices for SWIRE 
required low Galactic cirrus emission, large contiguous areas 
and low foreground contamination such as bright stars and local 
galaxies. The observing strategy for the SWIRE survey involved 
mapping of the fields such that each point is covered by four 30- 
second exposures with the IRAC instrument and 44 four-second 



1 ELAIS Nl, ELAIS N2, ELAIS SI, Lockman, XMM-LSS, CDFS 

2 IRAC: 3.6yt/m, 4.5yt/m, 5.8 yum, 8.0yum. MIPS: 24yt/m, 70yumand 
160yum. 



exposures with the MIPS instrument. The SWIRE observations 
were taken throughout 2004 with over 800 hours of observations. 

After initial analysis of the SWIRE fields we opt to remove 
the ELAIS S 1 field from comparison, due to its lower coverage 
compared to the remaining five SWIRE fields. By doing this, 
we manage to cut at a deeper flux density without jeopardising 
the signal to noise, whilst still allowing our sample to be large 
enough to mitigate cosmic variance. 

2.3. Reduction 

Data reduction was performed using MOPEX (Makovoz & Khan 
2005) as described by Seymour et al. (2007); basic calibrated 
data are drizzled to reach an increased sampling of 1725 per pixel 
from the intrinsic pixel scale for MIPS 24yum frames of 2! 5 per 
pixel^. 

3. Source Extraction and Photometry 

3.1. Creation of Variance IVIaps 

Fully reduced data and corresponding coverage maps are used 
to produce accurate variance maps: the pixel flux densities are 
binned and plotted according to their coverage, Gaussian profiles 
are fit to the lower part of the resultant histograms (ignoring the 
positive skew due to sources), then finally RMS noise values 
are extracted for each coverage value. Plotting of the root mean 
square (RMS) noise values v^'. coverage (and using cr^ oc I /"W, 
where cr is RMS noise and is the weight, equivalent to 
coverage) allows the coverage maps to be converted directly to 
absolute variance maps. 

3.2. Detection and Extraction 

Source detection and photometry are undertaken using the 
SExtractor software (Bertin & Arnouts 1996). Parameters are 
set as to minimise false detections and to optimise the separate 
detection of sources in close proximity. The parameters chosen 
are included in Appendix A. The flux density from extended 
sources^ is taken to be the Kron flux density, with the aperture 
correction to the flux density being calculated from the effective 
Kron radius. Elsewhere, the resultant flux density from point- 
like^ sources is taken to be that calculated within a 5V25 aperture 
radius. Aperture corrections are calculated using the Spitzer 
MIPS data handbook (see the MIPS data handbook^, p.37). 

3.3. Flux Density and Region Cut 

Due to the dithering pattern adopted during the HzRG field 
observations for which the central region of each frame has a 
higher coverage, we restrain the present analysis to the central 
1.75' radius circular fields centred on the HzRGs, this rejects 
areas in the frame with lower coverage. 

To determine a suitable flux density cut to our data, we first 
calculate the 5V25 radius aperture errors for each frame. This 
is done by placing random apertures upon the frame (within 
1.75' radii circles), binning of source flux densities, and plotting 
the flux density distribution. We again perform a Gaussian fit 
to the lower part of the distribution, so as to avoid the positive 

^ note that the SWIRE data is rescaled to V!2 per pixel 

^ area > 100 pixels and stellarity < 0.8 

^ area < 100 pixels or stellarity > 0.8 

^ http://ssc.spitzer.caltech.edU/mips/dh/mipsdatahandbook3.3. 1 .pdf 
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Table 1. SWIRE fields. In Col. 1 are the field names, in Col. 2 are 
the areas analysed, in Col. 3 are the RMS values calculated for 5'.' 25 
radii apertures (not aperture corrected). Column 4 states the average 
number of sources per 1.75' radius circular cell. The total number of 
1.75' radius circular fields analysed in SWIRE is 9981 (area 26.7 deg^), 
totalling 72295 sources. 



Field Name 


Area 


RMS 


N 

^ ^ sources 




(deg2) 


(A/Jy) 


(1.75' rad) 


(1) 


(2) 


(3) 


(4) 


ELAIS Nl 


6.73 


28.0 


7.16 


ELAIS N2 


4.60 


28.3 


7.83 


CDFS 


7.69 


27.9 


7.16 


Lockman 


7.81 


29.6 


7.27 


XMM-LSS 


7.40 


34.3 


6.03 



skew due to flux from sources, and calculate the Gaussian width; 
this width gives us the RMS noise for each frame (RMS values 
are included in Tables 1 and 2). Our RMS values calculated for 
the SWIRE fields are smaller than those of Shupe et al. (2008) 
by a factor of ~ 1.35; however, our results are consistent with 
their signal-to-noise measurements. For comparison, the average 
RMS noise for the HzRG fields is 23.0 yuJy, while for the SWIRE 
fields the average RMS noise is 29.6 /iJy. We adopt a flux density 
cut of f2A^m = 0.3 mJy, corresponding to a minimum S /N - 5 
for the fields with the highest RMS, after aperture correction 
(both for our HzRGs and for SWIRE fields, this is limited by 
the SWIRE XMM-LSS fields). As wiU be shown in Sec. 4.5, 
the HzRG RMS values are found to be independent of 24 jim 
source densities, which is indicative of a conservative cut free of 
questionable sources. 

3.4. Densities in 1.75' circular cells 

Contiguous 1.75' circular cells are placed intermittently (without 
overlap) within the five SWIRE fields, totalling some 9981 cells. 
The 24 yum source densities are then calculated in each of these 
cells, with an imposed flux density cut of /24//m = 0.3 mJy. 
Source densities in our HzRG fields are calculated by centreing 
the 1.75' radius apertures on the radio galaxy locations. When 
undertaking a sources-in-cell count for our HzRG fields, the 
HzRG itself is removed from the statistics. Uncertainties are 
calculated from small number statistics, following the method 
of Gehrels(1986). 

4. Analysis 

4.1. Comparison with SWIRE data 

Figure 1 shows the source density distribution of the five SWIRE 
fields. The grey histogram represents the percentage of cells 
containing a given number of sources. The dashed line shows 
the Gaussian fit to the lower portion of this distribution, and 
exemplifies the positive skew in the SWIRE data due to clus- 
tering. Our values for the Gaussian fit are (N) = 5.28 with 
(Tn = 2.52. Superimposed upon this plot is a histogram of the 
HzRG field densities, whereby the histogram height represents 
the number of HzRG fields for a given number of sources. Along 
with the HzRG histogram the HzRG field names are included. 
We also overplot our criterion for protocluster candidacy of 
NhzRG > {N} + 3crN = 12.84, a cut which only 6% of SWIRE 
cells exceed. This comparison of our HzRG field data with 
the SWIRE fields shows such an overdensity in a significant 
minority (20 [32%]) of HzRG fields. The mean density of the 



3 

HzRG fields is 11.6 ± 6.3, translating to a mean overdensity 
(S = of 2.2 ± 1.2 over the reference SWIRE data, with 

SWIRE 

a Kolmolgorov-Smirnov probability that the two datasets are 
drawn from the same distribution of 1.2 x 10"^^. Individual field 
source densities are included in Table 2 



4.2. Comparison with IRAC 

G12 study the environments ofl.2<z<3.0 HzRGs and com- 
pare the densities of sources with IRAC [3.6] -[4.5] > -0.1 (AB) 
(within a certain flux density cut) found within one arcminute 
of the HzRG. The imposed colour cut selects against galaxies 
below z ~ 1.2, but does not distintinguish between passive 
and active galaxies (Papovich 2008). Comparing the main result 
of this paper, a significant fraction of our sources confirm the 
protocluster candidacies proposed in G12. The combined result 
of this comparison is shown in Fig. 2, where the 24yum source 
number in 1.75' radius circular cells is plotted against radio 
galaxy redshift. The dotted line, dark shaded region and light 
shaded region represent the SWIRE mean density, (N), Icr^ and 
Sctn deviations, all extracted from the Gaussian fit described 
above. Plotted star symbols represent those sources with an 
overdensity in our work of Smips ^ (N) + Scr^ and in G12 with 
an overdensity of Sirac ^ (N) + 2(7^. Plotted circles represent 
those remaining sources present in both catalogs, and the plotted 
plus symbols represent the fields which are not analysed in 
G12: 3C65 (z=1.176) and LBDS53W059 (z=1.432). Notes on 
protocluster candidates are included in Sec. 4.6 and shown in 
Fig. 7. 



Figure 3 shows the source field densities as measured in 
IRAC by G12 versus the 24 jim source densities measured here. 
A strong correlation between the two values is found, despite the 
wavebands probing diflferent populations. Of particular interest 
for further study may be those environments which are found 
to be significantly dense in one waveband but unremarkable or 
even underdense in the other, as these field abnormalities may 
enlighten us with regard to structure formation. The Spearman 
Rank correlation is p = 0.577, implying a statistically significant 
correlation (e.g. there is a 99.992% probability that the two 
parameters are correlated). 



As a check of the quality of our data, the RMS noise is plotted 
against the 24yum source density (Fig. 4). We expect there to 
be no obvious correlations between these values, and indeed no 
correlation is found. The Spearman Rank Correlation coeflftcient 
for these data is p = 0.08, implying only a 48.6% probability of 
a correlation between the values. 



4.3. Evolution with Radio Luminosity 

Plotting of 24 yum source density versus HzRG 500MHz lumi- 
nosity (Fig. 5) shows no significant correlation between the two 
parameters. Though no work in the literature has had a large 
enough sample to statically prove or disprove any correlation, 
suggestions that high luminosity radio galaxies lie within a 
more dense medium have been put forward (Venemans et al. 
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Table 2. HzRG field Densities. HzRG names, redshifts, RMS noise values (for 5'/. 25 radius apertures, not aperture corrected) and calculated 
background values in Cols. 1-4. Exposure time is stated in Col. 5. 500 MHz radio luminosities quoted in Col. 6 are taken from De Breuck et al. 
(2010). 1.75' radius cells source counts and each field overdensity (with respect to the calculated SWIRE 1.75' radius fields Gaussian peak, as 
per Fig. 1) along with their associated small-number statistical error in Cols. 7 and 8. Column 9 shows the field overdensities as a function of cr^, 
corresponding to the probability of finding that given number of sources in a blank, non targeted field. Fields classified as protocluster candidates 
(20) are in bold font. 



HzRG 


Redshift 


RMS Noise 


Background 


Exposure 


loglo(L500MHz) 


# Sources 


6±Eyy 


<5(o-) 






L"Jy] 


[mJy] 


[s] 


[WHz-1] 


([1.75' radius]) 


[NhzRg/Nswire] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 



3C356.0 

MRC0037-258 

3C368.0 

6C0058+495 

3C65 

3C266 

MRC021 1-256 
MRC0114-211 
LBDS53w069 
7C1756+6520 

7C1751-h6809 
LBDS53w091 
3C470 

MRC2224-273 

6C0132-h330 

MRC1017-220 

3C239 
3C294 

7C1805+6332 

6CE0820-h3642 

6CE0905-h3955 

MRC0324-228 

MRC0350-279 

MRC0152-209 

MRC2048-272 

PKS1138-262 

4C40.36 

TXS021 1-122 

USS1707-hl05 

USS1410-001 

LBDS53w002 

6C0930-h389 

MRC0406-244 

4C23.56 

MRC2104-242 

WNJl 115-^5016 

PKS0529-549 

MRC2025-218 

USS2202+128 

MG1019-h0534 

4C24.28 

4C28.58 

USS0943-242 

WNJ0747+3654 

B3J2330-h3927 

MRC0316-257 

MRC025 1-273 

WNJl 123-^3141 

6C1232-h39 

TNJ0205-h2242 

TNJ0121-hl320 

TXJ1908-h7220 

USS1243-h036 

WNJ191 1-^6342 

MG2144-hl928 

6C0032-h412 

4C60.07 

4C41.17 

TNJ2007-1316 

TNJ1338-1942 

8C1435-h635 



1.079 

1.100 

1.132 

1.173 

1.176 

1.275 

1.300 

1.410 

1.432 

1.480 

1.540 

1.552 

1.653 

1.679 

1.710 

1.768 

1.781 

1.786 

1.840 

1.860 

1.883 

1.894 

1.900 

1.920 

2.060 

2.156 

2.265 

2.340 

2.349 

2.363 

2.393 

2.395 

2.427 

2.483 

2.491 

2.540 

2.575 

2.630 

2.706 

2.765 

2.879 

2.891 

2.923 

2.992 

3.086 

3.130 

3.160 

3.217 

3.220 

3.506 

3.516 

3.530 

3.570 

3.590 

3.592 

3.670 

3.788 

3.792 

3.840 

4.110 

4.250 
A A^'l 



19.5 
22.1 
27.6 
24.1 
24.0 
17.3 
20.6 
18.8 
15.8 
18.7 
18.0 
12.6 
25.4 
27.6 
28.3 
23.0 
26.7 
25.2 
18.5 
23.3 
17.5 
16.2 
27.0 
20.6 
23.1 
9.2 
22.3 
25.2 
17.7 
33.7 
22.8 
23.5 
28.7 
32.3 
28.4 
19.2 
30.5 
25.9 
21.4 
22.3 
20.4 
22.3 
17.3 
31.5 
19.4 
23.1 
15.9 
25.9 
19.2 
26.5 
28.3 
26.0 
29.4 
19.4 
25.6 
21.4 
31.4 
23.1 
18.7 
32.0 
21.0 

0/1 Q 



0.76 
0.75 
0.73 
0.85 
0.89 
1.07 
0.79 
0.79 
0.73 
0.78 
0.78 
0.73 
0.83 
0.75 
1.18 
0.72 
1.07 
0.49 
0.78 
0.78 
0.78 
0.79 
0.78 
0.79 
0.76 
1.08 
0.77 
0.53 
0.79 
1.24 
0.66 
1.39 
0.76 
0.77 
0.78 
1.47 
0.79 
0.76 
0.98 
1.44 
0.76 
0.76 
0.66 
0.81 
0.69 
0.79 
0.79 
1.45 
1.10 
0.77 
0.78 
0.48 
1.45 
0.77 
0.85 
0.77 
0.77 
0.78 
0.76 
0.78 
0.76 

1 



450 
900 
450 
450 
900 
900 
900 
900 
2100 
450 
450 
2100 
450 
900 
900 
900 
900 
900 
450 
900 
900 
900 
450 
360 
900 
9300 
450 
900 
900 
900 
450 
900 
450 
450 
900 
900 
450 
900 
900 
900 
900 
900 
450 
900 
450 
900 
900 
900 
900 
900 
900 
450 
900 
450 
900 
450 
450 
900 
900 
900 
450 



28.35 
27.72 
28.52 
27.33 
28.63 
28.54 
27.78 
28.66 
26.30 
27.40 
27.46 
27.04 
28.79 
27.52 
27.64 
27.94 
29.00 
28.96 
27.78 
28.28 
28.17 
28.49 
28.25 
28.20 
28.72 
29.07 
28.79 
28.48 
28.63 
28.41 
27.78 
28.41 
29.03 
28.93 
28.84 
27.82 
29.16 
28.74 
28.54 
28.57 
29.05 
28.91 
28.62 
28.14 
28.33 
28.95 
28.54 
28.51 
28.93 
28.46 
28.49 
29.12 
29.23 
28.14 
29.08 
28.75 
29.20 
29.18 
29.13 
28.71 
29.40 



7 
11 
17 
15 
15 
12 
16 
20 
15 
9 
10 
16 
11 
9 
15 
10 
9 
21 
7 
8 
7 
13 
12 
11 
16 
13 
8 

10 
9 
14 
7 
9 
17 
12 
8 

18 
27 
13 
10 
10 
10 
9 
15 
11 
16 
4 
4 
8 
6 
5 

12 
12 
7 
12 
11 
9 
11 
9 
10 
9 

1 Q 



1.5 
1.3 
2.1 
3.2 
2.8+1 
2.8+; 
2.3+1 
3.0 
3.8 
2.8+; 

1.7+; 
1.9+; 

3.0+ 

2.1+; 
1.7+; 
2.8+; 
1.9+; 

1.7 
4.0 
1.3 
1.5 
1.3 
2.5- 
2.3 
2.1 
3.0 
2.5- 
1.5 
1.9 
1.7 
2.6 
1.3 
1.7 
3.2 
23- 
1.5 
3.4 
5.1 
2.5 

1.9+; 
1.9+; 
1.9+; 
1.7+; 
2.8+; 

21 
3.0 
0.8 
0.8 
1.5 
1.1 
0.9 
23- 
23- 
1.3 
23- 
21 
1.7 
21 
1.7 
1.9 
1-7-0 
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Fig. 1. A comparison of HzRG source counts with the reference SWIRE fields. The distribution of the SWIRE 24 fjm source counts in ~ 10, 000 
cirular cells of 1.75' radius, representing five SWIRE fields, is plotted as a grey histogram. A Gaussian fit to the lower part of this distribution is 
overlayed as a dashed line. A 3cr above the SWIRE mean is indicated as a dot-dash line. A histogram of HzRG field densities is overlayed along 
with each HzRG name. The corresponding source number density is indicated on the right axis. 



2007; Miley & De Breuck 2008; Falder et al. 2010). A possible 
explanation for this is that higher luminosity radio galaxies lie 
in a medium with a higher than average density, thus being more 
conducive to cluster formation. Our work, alongside that of G12 
are the first capable of making such a statistically significant 
comparison and we find no such correlation between source 
density and HzRG 500 MHz luminosity. The Spearman Rank 
Correlation coefl&cient recovered for these data is p = 0.04, 
implying only a 25% probability that the two parameters are 
correlated. We note that there is a small redshift-radio power 
degeneracy in our sample (especially at z > 3, see Fig. 5). 
Given that redshift and density are clearly correlated (see Fig 2) 
this may imprint a correlation between source density and radio 
power. We take this into account by using the Spearman Partial 
Rank correlation coefl^icient (Macklin 1982) which reduces the 
correlation coeflftcient to p = 0.01. We conclude that the 24yum 
source density does not depend on radio power. 



4.4. K-correction and Population Type 

Our sample covers a wide range in redshifts, and the observed 
24yum band therefore covers significantly diflferent restframe 
wavelengths. To examine these K-correction eflTects, we map 
the Spectral Energy Distributions (SEDs) of a starburst galaxy 
(Siebenmorgen & Kriigel 2007) as well as face- and edge- 
on quasars (Pier & Krolik 1992) through the MIPS 24yum 
transmission window for redshifts 1 < z < 5.2 (Fig. 6). 

It is interesting to note a number of things; (1) the relative 
flatness of the quasar flux density over the entire redshift range, 
due to hot dust entering the MIPS transmission window; and 
(2) the bump in the starburst SED at z ~ 1.5, due to the 
PAH emission lines (dominated by 7.7 yum) entering the MIPS 
transmission window. Note that the smoothness of the starburst 
line is due to the broad MIPS 24yum transmission window 
covering some 9 jim between -20-29 yum, thus smoothing out 
the more narrow PAH features. Also note that even for such 
a high SFR starburst, galaxies fall below the flux density cut 
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Fig. 2. Source counts around each HzRG versus HzRG redshift. The 
dotted Une indicates the SWIRE mean source count for a 1.75' radius 
cell, the dark shaded region denotes +lcr and the light shaded region 
denotes ±3cr on the SWIRE data. Starred symbols denote those sources 
which are considered candidate clusters in G12 (8) while the plus 
symboled sources (2) are those which are not studied in G12. Note that 
z > 3 sources are not deep enough to clarify protocluster candidacy 
in G12. The uncertainties on the source counts are derived from small 
number statistics. 
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Fig. 3. Comparison of MIPS field densities with IRAC [3.6] - [4.5] > 
-0.1 field densities from G12 (1.2 < z < 3.0). The horizontal and 
vertical dotted line correspond to the mean SWIRE 24yL/m and IRAC 
colour- selected values from SWIRE, respectively. A relatively strong 
correlation between source densities is found. 




10 



15 



20 25 
RMS [uJy] 



30 



35 



Fig. 4. HzRG field density versus RMS noise. No obvious correlation 
between the two values is found. The error bar plotted in the upper left 
in included to give an indication of a typical uncertainty in the RMS 
values. SWIRE field values are plotted as larger circles where yellow, 
blue, red, green and cyan correspond to CDFS, ENl, EN2, XMM-LSS 
and Lockman fields respectively. 



beyond redshift three. Finally, the orientation of quasars will 
have a large impact on their detectability at lower luminosities. 

By comparing a high SFR starburst to a moderately luminous 
quasar, as shown in Fig. 6, we can make the assumption^ that 
above a redshift of z ~ 3 we are probing almost exclusively 
the quasar population. On the other hand, below z ~ 3, we 
are simultaneously probing the quasar and starburst populations, 
with starbursts become ever-more important at redshifts two and 
below. We are likely to see higher densities in fields at lower 
redshifts. This is indeed what we observe; as shown in Fig. 2 
the highest density fields (S(cr) > 4.5) almost exclusively lie 
below a redshift of three. Reddy et al. (2006, see Fig. 2) compare 
spectroscopically confirmed 24yum sources in the GOODS- 
North field between redshifts 1 < z < 2.6. For a flux density 
cut of f24j2m > 0.3 mJy, they find only Sub-lMillimetre Galaxies 
(S]V[Gs) and X-ray sources in the field, consistent with probing 
the starburst galaxy and AGN populations. 

The minimum AGN luminosities which we can probe with 
our imposed flux density cut are 3.4 x 10"^^ erg s"^ at z = 1 to 
1.9 X 10^^ erg s"^ at z = 5.2. Compare this with the minimum 
SFRs which can equivalently be probed; SFR= lOOMoyr"^ at 
z = 1 to a SFR -lOOOM© yr"^ at z = 3. These values have been 
calculated using the relevant formulae from Rieke et al. (2009). 

Since we do not know the physical nature of every object 
in each frame, let alone the redshift of each of these sources, 
it is not possible to perform K-corrections. For K-corrections 
to be calculated we will need to scrutinise the complimentary 
IRAC data for these fields, as well as gain data in the Far-IR. 
This process is currently being undertaken and will be the focus 
of future works. However, it is beyond the scope of this paper. 



if we assume no evolution of the luminosity functions. 
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Fig. 5. HzRG 24 yum field densities versus HzRG 500 MHz luminosities 
(taken from De Breuck et al. 2010). The HzRGs with z > 3.0 are plotted 
as open symbols. 



4.5. Evolution with Redshlft 

Figure 2 also shows the HzRG field densities versus HzRG red- 
shift. If 6C0140+326 (z=4.413) is a lensed system (as suggested 
in Lacy 1999) and as such possibly contaminated by artificially 
brightened sources, then this field should be removed from our 
sample. If we do remove this field, then we begin to see a general 
trend towards smaller-overdensities at early epochs. This could 
be interpreted in two ways: (1) the cluster luminosity function is 
truly evolving over time; or (2) the imposed flux density cut is 
too high to probe comparable populations at higher redshifts, as 
discussed above. 

Since we are probing diflferent populations at lower (z = 
1-2) and higher (z = 3 - 5) redshifts, with a transitory region 
between them, higher source densities at later epochs is the 
behavior that one would expect. However, since we also expect 
structure formation to be evolutionary, and not instantaneous, 
observations of higher overdensities at lower redshifts cannot 
be assigned to selection bias alone. The mid-IR luminosity 
function has been pushed out to z ~ 2.5 through stacking 
analysis (Caputi et al. 2007; Rodighiero et al. 2010) though 
is essentially unconstrained at z > 2.5. We do not expect a 
significant increase in number density at higher redshifts as 
this would be inconsistent with the star-formation rate density- 
evolutionary plot (Lilly et al. 1996; Madau et al. 1996); the 
high-end of the mid-IR luminosity function must drop oflT at 
high redshift, though constraints are too poor to infer population 
statistics. 
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Fig. 6. Observed MIPS 24/im flux density versus redshift. Solid 
black line represents a starburst galaxy with a star formation rate of 
1000 M© yr~^ . The dashed green line and dot-dash blue line represent a 
face-on and an edge-on quasar with L24^m = 10^^ erg s~^ . The horizontal 
line represents the 0.3 mJy flux density cut. 



4.6. Notes on Individual Sources 

For the sake of a well defined sample, we only select lAjim 
targets with an A^hzRG > (Nswire) + 3 cr as cluster candidates. 
This accounts for 20 (32%) of our sample. For all sources that 
meet this criterion, cutouts are included in Fig. 7. Fields are 
sorted by increasing redshift. Quoted sigma values correspond 
to the probability of finding protocluster candidates compared to 
a blank, non-targeted field. 



4.6.1. 6C0058+495 (z=1.173) 

We find a 24 jim source overdensity of 4.7(T, with the sources 
in this field having a clumpy structure. This field is our lowest 
redshift protocluster candidate. 



4.6.2. 3C65 (z=1.176) 

This field was first studied by Best (2000), who found a ^-band 
overdensity. Our study confirms an overdensity in 3C65 at a 3.9(T 
level. 



4.6.3. 3C266 (z=1.275) 

This field has also been studied in G12, but its source overdensity 
was not significant enough to declare protocluster candidacy. 
This work however finds a 24yum source overdensity of 3.9(T, 
making it a firm candidate. 



which concentrates only on the statistical properties of HzRG 
fields. 



4.6.4. MRC01 14-211 (z=1.41) 

With an overdensity of IRAC sources significant at a 6cr level 
compared to the field, MRCOl 14-211 is the densest field and 
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most promising cluster candidate in G12. We also find an 
overdensity of 24 fim sources at the 4.3 cr significance. 

4.6.5. LBDS53W069 (z=1.432) 

This is formally not a radio galaxy according to the radio lu- 
minosity cutoff' defined by Seymour et al. (2007). Nevertheless, 
our work shows a significant 24yum source overdensity at S.Scr; 
observationally this field is very dense with no preferred spatial 
segregation. This field was not included in G12's [3.6] - [4.5] > 
-0.1 analysis, but it has been calculated that this field would 
represent a 4.8(T overdensity. 

4.6.6. 7C1 756+6520 (z=1.48) 

This field is a known protocluster, first identified in Galametz 
et al. (2009) and spectroscopically confirmed in Galametz et al. 
(2010). This work confirms that 7C1756+6520 is also a proto- 
cluster candidate at 24 yum, with a source overdensity of 3.9(T. 

4.6.7. 3C470 (z=1.653) 

This field was first studied in G12, with a source overdensity of 
2.9(T. Our work quotes a 24 jim source overdensity of 43 cr. 

4.6.8. MRC1017-220 (z=1.768) 

This field was first studied by Cimatti et al. (2000), who found 
an overdensity of EROs. Studies of near-IR selected galaxies 
(Galametz et al. 2010) and IRAC-selected galaxies (G12), how- 
ever, did not find a significant overdensity. Our work finds a 
24 yum source overdensity at the 3.9cr level. 

4.6.9. 7C1 805+6332 (z=1.84) 

This field was first studied in Galametz (2010), where no near- 
IR overdensity was reported. This field was also studied in G12, 
which again did not find a significant overdensity. This work 
finds a 24/um source overdensity at 6.2(T, making it a very 
strong protocluster candidate. Observationally, sources lie in a 
filamentary structure in the NE-SW direction. 

4.7. MRC0350-279 (z=1 .9) 

Field studied in G12, with an IRAC source overdensity of 
~ 5.2(T derived. Here, the 24yum source overdensity is 3.1(T. 
Observationally the field is very much clustered towards the East 
of the HzRG. 

4.7.1. PKS1 138-262 (z=2.156) 

PKSl 138-262 is a well studied protocluster environment with 
overdensities of Lya and X-ray sources (Pentericci et al. 2000, 
2002), as well as Ua emitters and EROs (Kurk et al. 2004a,b). 
This work finds an overdensity of 24 jim sources at 4.3(T. 

4.7.2. 4C40.36 (z=2.265) 

This field was studied in G12, who found an underdensity of red 
IRAC sources. Our work quotes a 24 yum source overdensity of 
3.1(T. Sources appear to lie in a filamentary E-W structure. 



4.7.3. LBDS53W002 (z=2.393) 

This field was studied in G12, but no significant overdensity 
was found. The 24yum source overdensity quoted here is 3.5cr. 
The structure appears to have some form of preferred direction, 
though this remains ambiguous without a larger field of view. 

4.7.4. 4C23.56 (z=2.483) 

There are many works covering this field, beginning with Knopp 
& Chambers (1997) that quote a "marginally significant excess 
of objects" in the ^-band. Other works have found overdensities 
in the near-IR (Kajisawa et al. 2006), of Ha emitters (Tanaka 
et al. 2011) and mid-IR (G12). This work finds a 24yum source 
overdensity at a 4.1cr level. 

4.7.5. PKS0529-549 (z=2.575) 

This field was studied in G12 and its environment was found 
to be underdense. Here, however, a significant 24yum source 
overdensity of S.lcr is found. The sources appear to be very 
strongly clustered within the 1.75' radius frame, making this 
field a firm protocluster candidate at 24 fim. 

4.7.6. MRC2025-218(z=2.63) 

This is our highest density field, and has recently been studied 
by G12, quoting a 2.9cr overdensity. Our result shows a 24/jm 
source overdensity at 8.6(T, with 27 (/24^m > 0.3 mJy) sources in 
a 1.75' radius field. This exceeds the number of sources in every 
1.75' radius SWIRE sub-field from our analysis. 

4.7.7. USS2202+128 (z=2.706) 

This field was studied in G12, who found no overdensity. Here, 
the 24 yum source overdensity is quoted as 3.1(T making it a firm 
protocluster candidate. Clustering appears mostly in the North 
of the frame, centred on a bright object ~ 30" from the HzRG. 

4.7.8. WNJ0747+3654 (z=2.992) 

An 24 Jim overdensity of 3.9(T is found in this field, making 
WNJ0747+3654 a strong protocluster candidate. Upon visual 
inspection, local clustering of sources (within the 1.75' radius 
cell) appears to be taking place, though on larger scales (outside 
of the cell) no spatial segregation appears to be present. 

4.7.9. MRC0316-257 (z=3.13) 

This target is a known Lya and [QUI] protocluster (Venemans 
et al. 2005, 2007; Maschietto et al. 2008; Kuiper et al. 2010), 
which our work also confirms with a 24 yum source overdensity 
of4.3o-. 

4.7.10. 6C0140+326 (z=4.413) 

This is a known protocluster, having recently been studied by 
Kuiper et al. (2011), with a Lya overdensity of 9 ± 5 over a 
random field distribution. Our result confirms a 24 yum source 
overdensity in the vicinity of 6C0140-h326 at 5Acr. However, if 
it is lensed as suggested in Lacy 1999, then this field may be 
contaminated by artificially brightened sources. 
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5. Conclusions 

Mid-IR emission is a powerful tracer of both AGN (torus-dust) 
and starburst galaxies, making mid-IR observations ideal for 
finding galaxy clusters in the early Universe. Using targeted 
observations in the mid-IR, we have investigated the fields of 
63 HzRGs and find a statistically significant source overdensity 
over SWIRE reference fields data. We find a large minority of 
our targeted fields, 20, or 32%, to be overdense to at least a 
3(T significance, confirming 11 known protoclusters (or exist- 
ing protocluster candidates). We identify nine new protocluster 
candidates. Targeted mid-IR observations of HzRG fields shows 
to be a powerful technique for protocluster candidate selection. 

These overdensities are indicative of the presence of proto- 
clusters with active and star-forming galaxies at high redshifts. 
Our results indicate a redshift evolution of source density, but 
this may be due to an insensitivity to star-forming galaxies at 
high redshifts. Our dataset is compared with recent results of 
IRAC colour selected densities and correlations between the data 
are found. 

Further work regard these data will involve spectroscopic 
follow-ups for each of the cluster candidates in order to confirm 
members. Meanwhile Far-IR photometry currently being under- 
taken, in conjunction with the existing Spitzer data will allow 
a more complete knowledge of the population which make up 
these overdense fields. 



Appendix A: Photometry Parameters 

A.1. Basic SExtractor Parameters 

- CLEANing was turned off' 

- NPIXEL, the minimum number of adjacent pixels above the 
detection threshold to qualify a source, is set to 4. 

- DEBLEND_NTHRESH, the number of deblending sub- 
thresholds is set to 16. 

- DEBLEND_MINCONT, the minimum contrast parameter 
for deblending is set to 0.01. 

- PHOT_APERTURES, the fixed aperture diameter for pho- 
tometry is set depending on whether the SWIRE frames of 
the HzRG frames are being analysed. Since the re-sampling 
is diff'erent this value must be set to 8.75 or 8.4 respectively. 

- PHOT_AUTOPARAMS, The Kron-Factor and minimum 
Kron radius respectively are set to 2.55, 3.5. 

- THRESH.TYPE, is set to RELATIVE. This constrains 
DETECT.THRESH to be relative to the background. 

- DETECT.THRESH is set to 3. 

- BACK_SIZE, the background mesh size is set to 32 pixels. 

- BACK_FILTERSIZE, is set to 3. 

- BACKPHOTO.TYPE, is set to LOCAL. 

- SEEING_FWHM, the Stellar FWHM is set to 5.0 (arcsec- 
onds). 

- GAIN is set to 0.0 (e-/ADU) 

- WEIGHT.TYPE is set to MAP.VAR, since the weight map 
has been edited to be an absolute variance map. 

- WEIGHT.GAIN is set to N 
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Fig. 7. Cutouts of all HzRG 24 //m fields with a 3cr significance source overdensity. The dotted circle represents the 1.75' radius circular sub-field that is analysed 
further, and plotted in red are the source for which f24ijm ^ 0.3 mJy. All fields are scaled to the same 5.8' x 5.8' angular size, with North/East facing up/left. The 
HzRG is plotted as a blue star in each frame. Note that PKS1138-262 is much deeper GTO data. A colour counterpart is available in the online version of this article. 



